Localization and Practality in Inhomogeneous Quantum Walks with Self-Duality 
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We introduce and study a class of discrete-time quantum walks on a one-dimensional lattice. 
In contrast to the standard homogeneous quantum walks, coin operators are inhomogeneous and 
depend on their positions in this class of models. The models are shown to be self-dual with 
respect to the Fourier transform, which is analogous to the Aubry- Andre model describing the one- 
dimensional tight-binding model with a quasi-periodic potential. When the period of coin operators 
is incommensurate to the lattice spacing, we rigorously show that the limit distribution of the 
quantum walk is localized at the origin. We also numerically study the eigenvalues of the one- 
step time evolution operator and find the Hofstadter butterfiy spectrum which indicates the fractal 
nature of this class of quantum walks. 

PACS numbers: 05.60.Gg, 03.65.-w, 71.23.An, 02.90.-l-p 



I. INTRODUCTION 

A quantum walk (QW) ilhS], the quantum mechan- 
ical analogue of the classical random walk (RW), is a 
useful tool to promote various fields. For instance, the 
QW is an irnportant primitive in universal quantum 
computation ^ [sj and efficient quantum algorithms @- 
[l0| . The realization of topological phases [llj and quan- 
tum phase transitions in optical lattice systems by 
the QW have been discussed. The relation between 
the Landau-Zener transition and the QW was found in 
Ref. 13,]. The QW provides a new clue to more fun- 
damental problems such as the quantum foundations of 
time [14j] and the photosynthesis and efficient energy 
transfer in biomolecules [15), JJi] . There are several im- 
portant properties of the QW: (i) the inverted-bell like 
limit distribution, (ii) the quadratic speed-up of the vari- 
ance to the classical RW, and (iii) the localization due to 
the quantum coherence and interference. The physical, 
mathematical, and computational properties of the QW 
are summarized in Refs. (T7l-l20j. 

In this paper, we focus on the localization property of 
the QW and study a class of the discrete-time quantum 
walk (DTQW) on a one-dimensional lattice with spatially 
inhomogeneous coins. Throughout this paper, we define 
the localization such that the limit distribution of the 
DTQW divided by some power of the time variable has 
the probability density given by the Dirac delta function. 
The localization in the DTQW with homogeneous coins 
has so far been extensively studied. In the one dimen- 
sional DTQW, the localization was shown in the models 
with a three-dimensional coin [21.] . a four-dimensional 
coin [22], a two-dimensional coin with memory |23ll . a 
two-dimensional coin in a random environment [24| and 
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a time-dependent two-dimensional coin; the Fibonacci 
QW [1^. On the other hand, the nature of the local- 
ization in the two-dimensional DTQW has been stud- 
ied numerically [1^ [s^ and analytically [Sll, 111]. The 
DTQW on the Cayley tree with a multi-state coin was 
also studied js^] . Furthermore, the recurrence properties 
of the DTQW, i.e., the decay rate of the probability at the 
origin, were studied in Refs. (s^ - fssf . Recently, the local- 
ization property of the DTQW with the inhomogeneous 
two-dimensional coin has been analytically shown in the 
model in which the coin at the origin is different from the 
rest [1^ and in the model of the periodically inhomo- 
geneous coins [4l|. In this paper, we generalize the model 
investigated in Ref. jZlj to the incommensurate case and 
show the localization property of the generalized model. 

The generalized model defined in the next section is 
inspired by the Aubry- Andre model (42| , which describes 
the one-dimensional tight-binding model with an incom- 
mensurate potential. The corresponding Schrodinger 
equation is called the almost Mathieu equation (isj : 

ip{n + ip{n - 1) + V cos{2TTan)il]{n) = Eipin), (1) 

where ^{n) is the amplitude of the wavefunction at the 
nth site, a is the irrational number, and E is the eigenen- 
ergy. The model is often used to study the localization- 
delocalization transition [4^ ]. Aubry and Andre have 
shown that all the eigenstates are localized when V > 2 
while they are extended when V < 2. The model was 
shown to be self-dual with respect to the Fourier trans- 
form at V = 2 ^45]. The critical point V = 2 is closely 
related to the Azbel-Hofstadter problem, which is the 
eigenvalue problem of the Bloch electron on a square lat- 
tice in a uniform magnetic field [45l.|46}. The energy spec- 
trum shows a self-similar and fractal structure known as 
the Hofstadter butterfly [i^. The multifractal property 
of the spectrum and wavefunctions have been studied in 
Refs. 43, 3 and its Bethe ansatz solution was discussed 
in Refs. jU[5l|. In this paper, we show that our QW 
model is self-dual under the Aubry- Andre duality trans- 
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formation. We also study the eigenvalues of the one-step 
time evolution of our QW and find that the spectrum 
shows a self-similar and fractal structure in analogy to 
the Hofstadter butterfly. 

The paper is organized as follows. In Sec. we define 
the DTQW with spatially inhomogeneous coins and show 
that this model is self-dual under a certain transforma- 
tion, which is analogous to the Aubry- Andre duality. In 
Sec. mil we rigorously show that the limit distribution of 
our DTQW model is localized at the origin, which means 
that the distribution has the finite support for any time. 
In Sec. IIVI we numerically study the eigenvalues of the 
unitary operator for the one-step time evolution in our 
model and show the spectrum as a function of the inverse 
period (the Hofstadter butterfly diagram). Section Ivl is 
devoted to the summary and discussion. 



II. 



INHOMOGENEOUS QUANTUM WALK 
AND AUBRY-ANDRE DUALITY 



Let us define the one-dimensional DTQW (52| as fol- 
lows. First, we prepare the position and coin states de- 
noted a.s \n) ® \() = \n, () ■ Here, the position n is along 
the infinite one-dimensional lattice Z (the integers) and 
1^) is a normalized linear combination of \L) = (1,0)'^ 
and \R) — (0, 1)'^, where T is the transposition. Sec- 
ond, we describe the time evolution of the DTQW using 
the unitary operator U defined by two kinds of unitary 
operations C and W. The quantum coin fiip C G U{2) 
corresponds to the coin flip at each site in the RW. By 
the shift operator W, the position n changes depending 
on the coin state: 

W\n,L) := \n-l,L), 

W\n,R) := \n+l,R). (2) 

We can explicitly write the shift operator as 

ly = ^ (In - 1, L){n, L\ + \n + 1, R) {n, R.\) . (3) 

n 

Then, the one-step time evolution for the DTQW is ex- 
pressed as U = W{Id. dg) C), where Id. expresses the 
identity operator. Starting from the initial state |0, 0), 
we apply U to the state at every step. After the t time 
steps, we obtain the probability distribution at the posi- 
tion n as 

Pr(n;i)= ^ \{n,^\U'\0,<f,)\\ (4) 
i&{L,R} 

In the following, we introduce the DTQW with the 
spatially inhomogeneous coins. In this class of models, 
the quantum coin flip at each site can be different from 
the others. The coin operator C acting on the position 
and coin states is then given by 

C [(cos(27ran)|n, L) + sin(27rQ;n)|n, R)){n, L\ 

n 

-|-(cos(27rQ;7z)|n, R) — sin(27ran)|n, L)){n, R\] . (5) 



More explicitly, the operator C is written in the following 
form 



\ / I „ / cos(27rQ;7i) — sin(27ran) 
n)(n\ (Ki . )r, ( /o \ 
' ^ ' \ sm(27rQ;n) cos(27rQ!n) 



71 

:= [|n)(n| ®C'„ 



(6) 



where a is a real number (a G M) and corresponds to the 
inverse period of the coin operations. This inverse period 
can even be incommensurate to the underlying lattice. It 
reminds us of the Aubry- Andre model [l^j ■ We call this 
model the inhomogeneous QW throughout this paper. 
Note that, this model is a generalization of that in Ref. 
[4l| . As an example, in Fig. [T] we show the behavior of 
the probability distribution at the 1000th step obtained 
by numerical iteration. 




FIG. 1: Probability distribution in the inhomogeneous QW 
after the 1000 time steps. Here, we take the initial coin state 
(|L) + \R))/y/2 and the inverse period a — n/2, which is an 
irrational number. Note that, the probability outside the plot 
range is zero within the numerical accuracy. 

We shall now show that the inhomogeneous QW is 
always self-dual under a certain duality transformation. 
Let us introduce a duality transformation, which corre- 
sponds to the spatial Fourier transform [53|, 



n, L) = [sin(27rckmn)|m, L) + cos(27ramn)|m, R)] , 

rn 

n, R) ~ [cos(27ramn)|m, L) + sin(27ramn)|m, R)] .(7) 

m 

It should be noted that these states are not normalized. 
Let us now see the action of the shift operator W on 
these states. One can easily check the following relations 

W\n, L) = cos(27rQ!n)|ri, L) + sin(27ran)|n, R), 
W\n,R) ~ cos{2TTan)\n, R) — sin(27rari)|n, i). (8) 

It is now obvious that the shift operator W acts like a 
coin operator in the dual basis. Similarly, the action of 
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the operator C in the dual basis is given by 
C|n,L) = \n-l,L), 
C\n,R) = \n + l,R), 



(9) 



and it can be regarded as the shift operator in the dual 
basis. In this way, the duality transformation completely 
interchange the roles of W and C. Therefore, the in- 
homogeneous QW is self-dual under the transformation 



III. LIMIT LOCALIZED DISTRIBUTION 

In this section, we analytically show the localization in 
the inhomogeneous QW as the following theorem. 

Theorem 1. For any irrational a G R\Q, the limit dis- 
tribution of the inhomogeneous QW divided by any power 
of the time variable is localized at the origin: 

where Xt is the random variable for the position at the t 
step [5j^ , Ch. 2], 6 (> 0) is an arbitrary parameter, and 
"=4>" means convergence in distribution [55]. Here, the 
limit distribution I has the probability density function 
f{x) = d{x) (x € R), where S{-) is the Dirac delta func- 
tion. Note that, the limit distribution I is independent of 
the parameter 9. 



{t oo), 



(10) 



This theorem is the main result of this paper and 
means that the distribution of the inhomogeneous QW 
has a finite support for almost all a. We first show the 
following lemma to prove the above theorem in the case 
of specific rational a. 

Lemma 1. When a — (£ Q with relatively prime P 
(odd integer) andQ, the inhomogeneous QW is restricted 
to the finite interval [—Q,Q]. 

Proof. Let us express the state at the t step evolving from 
the initial state |0, 4>) as 



{wcY\o,^) = 



E 



(11) 



Then, we obtain the one-step time evolution of the coef- 
ficients 4'ti'n;^) as follows: 

V't+i(n - l;L) 

— cos{2T:an)'ipti'n-] L) — sin(27ran)'(/'t(n; R), 
V't+i(n + l;i?) 

= sin(27ran)-0t [n; L) + cos(27ran)-0t {n; R) . (12) 

Since cos(±27ramQ) = for any odd integer m, that is, 
the diagonal elements of the coin operator are zero at 
position ±mQ, we immediately arrive at 

V^t+i(±mQ -1;L) = ±{-1)'^ M±rnQ- R), 
i;t+i{±mQ + l;R) = ^{-1)"^ M±mQ; L). (13) 



The above equations mean that the quantum walker is 
reflected at the position n = ±mQ and cannot be trans- 
mitted across this point (see Fig. [2|). Since the initial 
position is localized at the origin (n — 0), the quantum 
walker undergoes reflection only at the sites ±Q. There- 
fore, we conclude that the distribution of the inhomoge- 
neous QW has the finite support [— Q, Q] for any time t. 

□ 



±mQ — 1 ±mQ 
site — I 1 — 




— I 1 




step 



FIG. 2: Time evolution of the inhomogeneous QW around 
the site zbmQ at which the quantum walker is reflected. The 
coin operator C and shift operator W act on the wavefunction 
\^n{t)) — {ipt{n, L),ipt{n, R))'^ at each site n. The wavefunc- 
tion around the reflection point ±mQ is indicated from the 
{tr — l)th step to the {tr -\- l)th step. 



From Lemma 1, we can conclude that our inhomoge- 
neous QW is finitely restricted when a — G Q. We 
now try to extend this result to the case of a e R \ Q and 
show Theorem [T] The following well-known theorem in 
number theory is useful. The proof of this theorem can 
be found in Ref. (56| . 

Theorem 2 (Dirichlet's Approximation Theorem [sgI . 
Theorem 185]). For any irrational number a G M \ Q, 
there is an infinity of the fractions a/b which satisfy 



< 



1 

62- 



(14) 



Theorem [2] implies the following corollary which will 
be used to prove Theorem 1. 

Corollary 1. For any irrational number a G R\Q, there 
is an infinity of the fractions with the relatively prime 
P ( odd integer) and Q which satisfy 



P 
4Q 



< 



4g2 



(15) 
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Proof. From Theorem [2l for any irrational number a' £ 
R \ Q, there is an infinity of the fractions P' /Q which 
satisfy 



P' 



< 



(16) 



P' can be expressed as P' = P x 2^ with an odd integer 
P and a non-negative integer £. It is noted that P and Q 
are relatively prime since P' /Q is an irreducible fraction. 
Therefore, dividing both sides of the above equation by 
we obtain 



Since a := a'/2^ 
desired result. 



P 
4Q 



< 



1 



< 



1 

4Q2- 



(17) 



must be irrational, we obtain the 

□ 



In the following, we prove Theorem [TJ 

Proof of Theorem]^ For any e > 0, there exist the rela- 
tively prime P (odd integer) and Q > 57 such that 



P 
4Q 



(18) 



where -1/(4(5^) < 6 < 1/(4(5^) according to Corollary[T] 
We calculate the diagonal elements of the coin operator 
at the ±Qth site as 



cos 



27r 



(±Q) 



= (-l)(^+i)/2sin(2^Q5). 

(19) 

Note that both the diagonal elements of C±q are equal. 



We obtain 



(-1) 



(P+l)/2 



sin(27rQ(5) 



< 



2Q 



< e. 



(20) 



Since e can be arbitrarily small, the diagonal element is 
zero in the limit. Then, along the same lines as Lemma[TJ 
the inhomogeneous QW is restricted to the finite interval 
[—Q,Q]. Since the time variable t and the parameter e 
are independent, the time variable t and the parameter 
Q are also independent. Therefore, we obtain the desired 
resuh (Hoi. □ 



Analogous to the above discussion, we show the lo- 
calization property of the DTQW with spatially random 
coins as follows. From Lemma [TJ we obtain the following 
remark, which was implicitly shown in Ref. [411. 

Remark 1. If the diagonal elements of the coin operator 
at a position on each side the origin are zero, then the 
DTQW is finitely restricted, i.e., the localization at the 
origin as in Eq. IJOj) . 

Therefore, the limit distribution of the DTQW with a 
randomly chosen quantum coin Cn,rsindom 

with uniform 

distribution in 17(2) at each position n (e Z) divided by 
any power of the time variable is also localized. This is 
because the diagonal elements of Cn, random must be zero 
at some positions. This result is relevant to Ref. [24 1. 



IV. FRACTAL PROPERTY 

In this section, we study the eigenvalues of the operator 
WC, which is the one-step time evolution in the inhomo- 
geneous QW and is defined in Eqs. ([31[5]). We show the 
general properties of the eigenvalues as well as the nu- 
merically obtained distribution of them, which suggests 
the self-similar and fractal structure in the spectrum in 
analogy to the Hofstadter butterfly problem. 

We first show the following theorem for the eigenvalues 
of the operators WC and CW. 

Theorem 3. All the eigenvalues of the operators WC 
and CW are identical including the degeneracy. 

Proof. Let be an eigenvector of the operator WC with 
a nonzero eigenvalue /i: 



WC\fi)^^,\ti). 



(21) 



Note that, all the eigenvalues of WC are nonzero and 
have modulus 1 since WC is a unitary operator. Multi- 
plying C from the left, one obtains 



CW{C\fi)) = ^^{C\^^)). 



(22) 



Since C is also a unitary operator, C|/i) is a nonzero 
eigenvector of CW with the eigenvalue ^. To complete 
the proof, we need to repeat the same argument for CW. 
Let I A) be an eigenvector of CW with an eigenvalue A: 



CW\\) = A|A). 
Multiplying W from the left, one obtains 
WC{W\\)) = \{W\X)). 



(23) 



(24) 



Since W is a unitary operator, W\X) is a nonzero eigen- 
vector of WC with the eigenvalue A. The identical de- 
generacy of WC and CW follows from the fact that the 
unitary operator C or preserves the orthogonality of 
the states with the same eigenvalue. □ 

Remark 2. Theorem \^ holds not only for the inhomo- 
geneous QW but also for the standard DTQW, in which 
the coin operator does not depend on the positions. 

According to Theorem [3J we have only to study the 
eigenvalues of CW. In the following, we focus on the case 
of the rational number a — & Q. In this case, from 
Lemma [TJ we can express CW as a 4(5 x 4Q irreducible 
matrix. Let us first derive a finite-dimensional matrix 
representation of CW. We can express the wavefunction 
at the t step evolving from the state |0, (p) by CW: 



(25) 



Along the same lines as the proof of Lemma [1] one 
obtains the one-step time evolution of the coefficients 
ipt{n, ^) as follows: 



ipt+iin;L) 
(pt+i{n;R) 



= Cn 



'fit (n - 



1;L) 
l;i?) 



(26) 
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where Cn is defined in Eq. ([6]). 

Since the initial condition is |0,(/)), one can clearly see 
that (pt{—Q] L) = (pt{Q\ R) — for any time t and hence 
the inhomogcncous QW by CW is also restricted to the 
finite interval [—Q,Q] like Lemma [TJ It is convenient to 
introduce the vector ipt defined by the coefficients ipt (n; ^) 
as 

0t = {ipt{-Q;R),(pti-Q + ^-^L),ipti-Q + 1;R), 

ifitiQ ~ 1; L),iptiQ - 1; i?), (^t(Q; L))'^. (27) 

Note that, ipt is the 4Q-dimensional vector. Using (pt, we 
can simply express Eq. (^5]) as 



ipt+i = CW(^t 



(28) 



Here, C and W are AQ x 4(5 matrices and their explicit 
forms are given by 



c = 



Q+l 







Q+2 







(29) 



W = 

where 
= 



/ f 

62 O ei 

62 O 












ei 



62 



ei 



O 



61 

/ 



(30) 



62 



(31) 



We note that empty entries in the matrices are zeros. One 
can easily confirm that C and W are unitary matrices. 

We shall now study the eigenvalues of CW. Since the 
matrix CW is unitary, the eigenvalues A(s C) have modu- 
lus 1, i.e., |A| = 1. Therefore, the argument of A faithfully 
represents the eigenvalues. Figure [3] shows the spectrum 
as a function of the inverse period a, obtained by ex- 
act diagonalization of CW. Although we can only study 
the case of rational a = ^ £ Q by numerical diago- 
nalization, we can understand the spectrum of CW with 
a G M \ Q as the irrational limit of a well-organized se- 
quence of rational numbers, whose existence is guaran- 
teed by Theorem [2l From Fig. [3l one can clearly see 
the self-similar and fractal nature in the eigenvalue dis- 
tribution. We also notice the following properties of the 
spectrum: 

(PI) All the eigenvalues at a are identical to those at 
1 — a. 

(P2) The eigenvalues come in complex conjugate pairs: 
for every eigenvalue A, there is an eigenvalue A*. 




FIG. 3: The self-similar and fractal structure in the inhomo- 
geneous QW. Arguments of the eigenvalues of CW (vertical 
axis) are plotted as a function of the parameter a = (hor- 
izontal axis) with Q < 50. Note that, P and Q must be 
relatively prime. 



(P3) The eigenvalues come in chiral pairs: for every 
eigenvalue A, there is an eigenvalue —A. 

(P4) All the eigenvalues are simple, i.e., nondegenerate. 

(P5) There are four eigenvalues A — ±1, ±z for any a — 

Let us explain the above properties in more detail. 
(PI) simply means that the spectrum in Fig. [3] is sym- 
metrical with respect to the line a = 1/2. This property 
can be shown as follows. Suppose that (f is the eigenvec- 
tor of CW at the inverse period a with the eigenvalue A. 
Then, according to Eq. (|26p , one obtains 



Xpi-Q;R) - (-1) — (^(-g + l;L), 



ip{n;L) 
Lp{n;R) 



= Cn{a) 



(p{n + 1; L) 
(p{n — 1; R) 



\p{Q-L) = {-l) — ^{Q-l-R), 



,inei-Q,Q)) 
(32) 



where C„ is denoted as C„(ck) to emphasize the argument 
a. We now consider the following transformation: 

ip{n; L) ^ (p{—n; R), (p{n; R) ^ (p{—n; L). (33) 

Under this transformation, Eq. ((32)) becomes 



Xvi-Q;R) - (-1) — (^(-Q + 1;L), 

Lp{n + 1; L) 
p{n~l;R) 



Xcp{Q;L) = {-l)^^iQ-l;R). 



.{n€{-Q,Q)) 
(34) 



The above equations are equivalent to CVJif — Xip at the 
inverse period 1 — a. Therefore, the sets of eigenvalues at 
a and 1 — a are identical. Furthermore, one can construct 
the eigenvectors of CW at 1 — a from those at a using the 
transformation ([55)1 . (P2) implies that the spectrum is 
symmetrical with respect to the line arg(A) = in Fig.|31 
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On can easily show (P2) by taking the complex conjugate 
of Eq. ([5^ and see that 0* is the eigenvector of CW with 
the eigenvalue A*. (P3) combined with (P2) gives rise to 
the symmetry with respect to the line arg(A) = n/2 or 
arg(A) — — 7r/2 in Fig. [3] (P3) is reminiscent of the chi- 
ral or particle-hole symmetry in the original Hofstadter 
problem, i.e., the eigenenergies of the Hamiltonian come 
in pairs ±E. This property holds not only for our inho- 
mogeneous QW but also for rather generic classes of the 
DTQW. To show (P3), we consider the local gauge trans- 
formation; (^(n;^) — )■ (— l)"(/3(n; ^). Under this trans- 
formation, Lp becomes the eigenvector of CW with the 
eigenvalue —A according to Eq. ([32)) . This proves (P3). 

Properties (P4) and (P5) requires a more elaborate 
analysis. To show (P4), we first show that tp{—Q; R) ^ 
if is an eigenvector of CW. Suppose for the sake of 
contradiction that (p{—Q;R) — 0. Using Eq. ([32)) . it 
follows that (f = 0. But this contradicts the fact that (f 
is the eigenvector of CW. Therefore, if{~Q\ R) ^ 0. It 
should be noted that 27m^ t:1/2 (i £ Z) for -Q < 
n < Q. We next suppose that CW has two eigenvectors 
0^^") and (^^^^ with the same eigenvalue A. Then, we define 
the linear combination of them as 

^ = ^(2) (_Q; i?)<^(l) - <p(l) (_Q; i?)<^(2) (35) 

and find 0'{—Q]R) — 0. The vector 0' is also an eigen- 
vector of CW and this contradicts the previous statement. 
Therefore, all the eigenvalues of CW are simple. Let us 
finally show the last property (P5). We first note that 
the product of the eigenvalues Afe [k — 1, 4Q) is 

4Q 

n^fe==-l (36) 

k=l 

since det CW ~ det(C) dct(W) = — 1. Here, we have used 
the fact that detC ~ 1 and detW = —1. This can be 
shown by an explicit calculation using Eqs. (j29|) and 
(150)) . Next, from (P2) and (P3), we see the eigenvalues 
Afe come in quadruplets Afe, A^, —Afe, — A^. if Afe is neither 
±1 nor ±i. Since |Afep — 1 and Eq. (|36)) . there exists k 
such that Afe — —1. The multiplicity of Afe — —1 should 
be 1 due to (P4) . Then, noting the fact that the number 
of eigenvalues 4Q is a multiple of 4, we can conclude that 
Afe ~ ±i should also be included as the eigenvalues of 
CW. This proves the statement (P5). Note that, the 
properties of the chiral symmetry (P3) and the existence 
of eigenvalues A = ±1 are relevant to Ref. [ll| discussed 
in the standard DTQW. 



analytically shown that the limit distribution of the in- 
homogeneous QW divided by any power of time is lo- 
calized at the origin when the inverse period of coins a 
is irrational number. We emphasize that our analytical 
solution for the localization is powerful since the round- 
off error is inevitable in numerical simulations. We have 
also numerically studied the eigenvalue distributions of 
the quantum walk operator WC in the inhomogeneous 
QW with a — G Q. The obtained spectrum shows 
the self-similar and fractal structure similar to the Hofs- 
tadter butterfly. 

In our model, the following problems remain open: 
First, it might be possible to experimentally realize the 
inhomogeneous QW in quantum optics as pointed out in 
Ref. [57|. Second, the recurrence property of the inho- 
mogeneous QW has not been analytically obtained. The 
example of the recurrence is numerically shown in Fig. |4| 
Finally, we have not yet obtained the limit distribution 
in the case of the rational a G Q except for a — 
with relative prime P (odd integer) and Q. According to 
Ref. ^58^, the limit distribution of the DTQW is partly 
related to the eigenvalues of the quantum walk opera- 
tor WC. From Fig. |3l the eigenvalues in the considering 
case seem to be continuous. Therefore, we propose the 
following conjecture as a general property of the DTQW. 

Conjecture 1. The distribution of the DTQW is not 
localized if the quantum walk operator WC only has con- 
tinuous spectra and does not have embedded eigenvalues. 

This conjecture is motivated by Ref. [H^ to show the 
property of the bound states from the spectrum of the 
discrete Schrodinger operator, which includes the almost 
Mathieu equation (HJ. Let us examine this conjecture for 
some cases in our inhomogeneous QW. It is easy to check 
that, in the case of a = P/2 (P G Z), WC has a contin- 
uous spectrum and localization does not occur. In this 
case, the coin state is unchanged up to the phase factor 

since the coin operator is always ^ q ^ '''^ ^ ^1 ) 

Therefore, the quantum walker does not undergo reflec- 
tion at any point. Furthermore, in the case of a = 1/6, 
Linden and Sharam showed that the standard deviation 
is proportional to the time variable t for a sufRcient large 
t [4l|. These facts partially support the above conjecture. 
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FIG. 4: Recurrence property of the inhomogeneous QW with 
the initial coin state (|L) + |i?))/\/2 and the inverse period 
a = 7r/2. The numericaUy obtained probabihties at the origin 
are shown up to 1000 time steps. Note that the probability 
at the origin must be zero when the step t is odd. 
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